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ABSTRACT 

One of the leading models for electron acceleration in solar flares is stochastic acceleration by weakly turbulent 
fast magnetosonic waves (“fast waves”). In this model, large-scale flows triggered by magnetic reconnection excite 
large-wavelength fast waves, and fast-wave energy then cascades from large wavelengths to small wavelengths. 
Electron acceleration by large- wavelength fast waves is weak, and so the model relies on the small-wavelength waves 
produced by the turbulent cascade. In order for the model to work, the energy cascade time for large-wavelength fast 
waves must be shorter than the time required for the waves to propagate out of the solar-flare acceleration region. To 
investigate the effects of wave escape, we solve the wave kinetic equation for fast waves in weak turbulence theory, 
supplemented with a homogeneous wave-loss term. We find that the amplitude of large- wavelength fast waves must 
exceed a minimum threshold in order for a significant fraction of the wave energy to cascade to small wavelengths 
before the waves leave the acceleration region. We evaluate this threshold as a function of the dominant wavelength 
of the fast waves that are initially excited by reconnection outflows. 
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1. INTRODUCTION 

Solar flares emit photons across a broad range of wavelengths, 
with ~99% of the radiated energy in the optical and UV contin- 
uum (Woods et al. 2004; Emslie et al. 2005). Although X-rays 
and gamma-rays account for only a small fraction of a flare’s 
total radiative luminosity, they provide important diagnostics 
on the populations of high-energy particles produced by flares. 
For example, flare X-ray spectra are typically thermal at pho- 
ton energies <20 keV and non-thermal (power laws or broken 
power laws) at energies >20 keV (Miller et al. 1997). Most 
of the X-ray emission is likely bremsstrahlung radiation re- 
sulting from Coulomb collisions between thermal protons in 
the chromosphere and high-energy electrons streaming down 
from the corona (Aschwanden 2002; Petrosian et al. 2002). 
Since bremsstrahlung X-ray photons of energy E\ are typically 
emitted by proton-electron collisions involving electrons with 
energies ~£'i, the observed X-ray spectra imply that the elec- 
trons typically have thermal distributions at energies <20 keV 
and non-thermal distributions at energies >20 keV (Miller et al. 
1 997). The intensity of the X-ray emission can be used to deduce 
the electron acceleration rate. For example, during the peak of a 
typical, large X-class flare, the observed X-ray emission implies 
that electrons are accelerated to energies >20 keV at a rate of 
~10 37 s" 1 (Miller et al. 1997). 

The origin of these high-energy electrons is a long-standing 
puzzle. Three main explanations have been advanced in the 
literature: shock acceleration, acceleration by large-scale, co- 
herent electric fields, and stochastic particle acceleration (SPA; 
Miller et al. 1997). In this paper, we focus on the SPA model, 
in which particles are accelerated by electromagnetic fluctua- 
tions associated with waves and/or turbulence. In this model, 
magnetic reconnection in the corona triggers rapid plasma out- 
flows from the reconnection site (Carmichael 1964; Hirayama 
1974; Kopp & Holzer 1976; Tsuneta et al. 1992; Tsuneta 
1996; Priest & Forbes 2000). Sunward flows encounter closed 
magnetic loops lower in the solar atmosphere, generating a 
disordered, turbulent flow above the loop tops. Initially, the 


velocity fluctuations in this disordered flow have a correlation 
length that is some fraction (perhaps ~1/10) of the loop size 
(~ 1 0 9 cm). This length scale is much larger than the proton in- 
ertial length i>a/O p = 2.3 x 10 2 cm- («/10 10 cm -3 ) -1 / 2 , where 
v A = Bq j sj An p is the Alfven speed, Bq is the strength of the 
background magnetic field, n is the proton number density, p is 
the mass density, and Q p is the proton cyclotron frequency. As 
a consequence, the fluctuations can be approximated as a super- 
position of magnetohydrodynamic (MHD) waves, i.e., Alfven 
waves, fast magnetosonic waves (“fast waves”), and slow mag- 
netosonic waves. 

Most SPA models of electron acceleration focus on inter- 
actions between electrons and fast waves through a type of 
resonant wave-particle interaction called transit-time damping 
(TTD; Miller et al. 1996; Schlickeiser & Miller 1998; Selkowitz 
& Blackman 2004; Yan & Lazarian 2004; Yan et al. 2008). 
In TTD, waves and particles satisfy the Landau resonance 
condition 

co/k\\=v h (1) 

where w is the wave frequency, k\\ is the component of the wave 
number k parallel to the mean magnetic field B o, and iiy is 
the component of the particle velocity along Bq. Equation (1) 
implies that particles “surf” on the wave phase fronts. In TTD, 
surfing particles are accelerated along Bq by the magnetic- 
mirror force arising from the perturbations in the magnetic field 
strength associated with the waves (Stix 1992). 

Another ingredient in the electron SPA models cited above 
is turbulence. After large-scale fast waves are excited by 
reconnection outflows, the waves interact with one another, 
causing wave energy to cascade from large wavelengths to 
smaller wavelengths. This energy cascade plays an important 
role in SPA models, because the small-wavelength fast waves 
produced by the turbulent cascade make the largest contribution 
to the electron acceleration rate via TTD (Miller et al. 1996). 

Our principal aim in this paper is to investigate how elec- 
tron acceleration in SPA models is affected by the escape 
of fast waves from the flare acceleration region. In low-/! 
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plasmas, where 


is the ratio of the thermal pressure p to the magnetic pressure, 
the fast-wave dispersion relation is 

cu = kvA- (3) 

Equation (3) implies that fast waves can propagate at speed Va 
in any direction. Fast waves can therefore cross magnetic field 
lines and leave the acceleration region. If the time required for 
waves to propagate out of the acceleration region is shorter than 
the energy cascade time at large wavelengths, then the large- 
wavelength waves that are excited by reconnection outflows 
will leave the acceleration region before their energy cascades 
to small wavelengths, sharply reducing the electron acceleration 
rate. This point was made previously by Miller et al. (1996). 
However, these authors did not attempt to calculate the degree 
to which the electron acceleration rate is reduced by wave escape 
as a function of flare parameters. In addition, their treatment of 
the fast-wave cascade was based on an estimate of the energy 
cascade time that was accurate only to within an unknown factor 
of order unity. We revisit the effects of wave escape using weak 
turbulence theory to calculate the energy cascade time from 
first principles. We describe our model of fast-wave turbulence 
and wave escape in Section 2, present our numerical results in 
Section 3, and discuss our conclusions in Section 4. 

2. LEAKY-BOX MODEL OF FAST- WAVE 
TURBULENCE IN FLARES 

We employ a “leaky-box” model to describe turbulence 
and wave escape in flares. We treat the flaring region as a 
homogeneous plasma with a uniform magnetic field Bq and 
define the fast-wave power spectrum F(k), abbreviated F*, to 
be twice the fast-wave energy per volume in Uspace per unit 
mass. The total fast-wave energy per unit mass is thus given by 

U t =^fF k d 3 k. (4) 


The last term on the right-hand side of Equation (5) provides 
a simple, hyperviscous model for the sink of fast-wave energy 
at large k resulting from TTD, which transfers energy from fast 
waves to electrons in flares, as described in the introduction. 
We have made this term oc k s purely for convenience, so that 
it “turns on” only at the largest wavenumbers in the numerical 
simulations that we present in the next section. The reason we 
include a factor of sin 2 0 in this term is discussed below. 

The form of the term (3E)./3f) tur b in Equation (5) depends 
on whether the turbulence is weak or strong. Fast waves at 
wavevector k are weakly turbulent when 

Sv k «; v A , (8) 

where Svk is the rms amplitude of the fast- wave fluctuation ve- 
locity at length scale k~ l (Kadomtsev 1965). When Equation (8) 
is satisfied, the linear fast- wave period P = 2 n/kvj^ is much 
shorter than the energy cascade timescale r c at wavevector k, 
and the fluctuations can be described as waves to a good approx- 
imation. If Equation (8) is not satisfied, and instead §Vk > i>a, 
then the turbulence is strong at wavevector k, and r c < P. 

Observationally, it is not clear whether solar-flare acceleration 
regions are in the weak-turbulence or strong-turbulence regimes, 
and in this paper we do not attempt to resolve the uncertainty 
in the wave amplitudes in solar flares. Instead, we assume 
that Equation (8) is at least marginally satisfied at all k (i.e., 
that 8vk ^ 0.3i>a at all k), so that weak-turbulence theory is 
approximately valid. We focus on this limit because wave escape 
is most important in this parameter regime. We will not model 
quantitatively the strong-turbulence case in Section 3. However, 
our results show that if the wave amplitudes are increased so 
that they approach the strong-turbulence regime, then the energy 
cascade proceeds sufficiently rapidly that most of the fast-wave 
energy dissipates before the waves can escape from the solar- 
flare acceleration region (see Figure 5). 

The wave kinetic equation for weakly turbulent fast waves 
in low-/l plasmas was derived by Chandran (2005, 2008). We 
set (3Ei/3r) tur b equal to the right-hand side of Equation (6) of 
Chandran (2005), with the Alfven-wave power spectrum A k in 
that equation set equal to zero: 


For simplicity, we assume reflectional symmetry, F(—k) = 
F(k). We take F k to evolve in time according to the equation 


3 t 




+ S k — sur(9)k s vFk, 
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where 6 is the angle between k and B 0 . The first term on the 
right-hand side of Equation (5) models wave escape from the 
flare acceleration region, which is characterized by the timescale 


_ h 

Tesc — ' 

va 


(6) 


where Lf is the size of the flare acceleration region. The term 
(3 Fk/dt\ m \, in Equation (5) gives the contribution to dFk/dt 
from wave-wave interactions, discussed below. The next term. 


S k 


4E 0 ( k\ 2 ( k 2 \ 

3 6XP \ k 2 )’ 


(7) 


is a source term representing fast-wave injection from recon- 
nection outflows. Here, Eq is the total wave energy injection 
rate per unit mass, and ko is the wavenumber at which Sk peaks. 


We have neglected Alfven waves for simplicity but expect that 
their inclusion would not significantly alter our conclusions 
about the effects of wave escape on fast-wave turbulence. 
Equation (9) is equivalent to the wave kinetic equation for weak 
acoustic turbulence, up to an overall coefficient (Zakharov & 
Sagdeev 1970; Zakharov et al. 1992). In both weak fast-wave 
turbulence and weak acoustic turbulence, waves with collinear 
wavevectors k , p, and q that satisfy the wavenumber resonance 
condition k = p + q and frequency matching condition k = 
p + q interact to produce a weak form of wave steepening, which 
transfers wave energy from small k to large k along radial lines 
in fc-space. However, fast-wave turbulence differs from acoustic 
turbulence in the following way. As sin 0 decreases — i.e., k 
and Bq become more closely aligned — fast waves become 
less compressive, the fast-wave cascade weakens, and the 
energy cascade time increases. This anisotropy is represented 
mathematically by the coefficient of sin 2 6 in Equation (9). It is 
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Figure 1. Steady-state fast-wave power spectrum at 9 = 50° in three numerical 
solutions in which koLf = 22. The thin solid line shows the scaling k 1 F k oc 
k~ 3 / 2 that arises in weak turbulence calculations without wave escape. 



Figure 2. Star symbols are the decay timescales T c computed from numerical 
integrations of Equation (5) for decaying fast- wave turbulence in the absence of 
wave escape (S k — > 0 and T esc -> oo), and the solid line is an analytic fit to the 
numerical results, given by Equation (12). 


also why we have made the hyperviscosity term in Equation (5) 
proportional to sin 2 6: this dependence prevents the dissipation 
term from truncating the fast-wave cascade at very small k when 
sin 6 is small. 


3. NUMERICAL RESULTS 

We numerically integrate Equation (5) forward in time 
until a steady state is reached. To evaluate the nonlinear 
term (d F k /dt) m b, we use the energy-conserving numerical 
algorithm employed by Chandran (2005), which is an extension 
of the numerical method developed by Leith & Kraichnan 
(1972) for statistical turbulence theories. We use a logarithmic 
wavenumber grid in both k± and k\\ (the components of k 
perpendicular and parallel to Bo), with k±j = ( 0.2ko)2 '/ 4 , 
i = 0, 1,2, ..., N - 1, fc|| 0 = 0, k\\j = (Q.2k () )2 (] - l)l \ 
j = 1, 2, . . . , N — 1, and N — 60. In all of our calculations, 
we choose the hyperviscosity coefficient v so that dissipation 
is negligible at small k (k < 10 2 ) but strong enough at large k 
(k > 10 3 ) to truncate the cascade. 

In Figure 1, we plot three steady-state solutions for F k at 
6 = 50°. For each of these solutions, k <} L( = 22. (As shown in 
the Appendix, the correlation length associated with an isotropic 
fast-wave spectrum localized at k = ko is L c = 2.2/ ko; thus, 
the choice koLf = 22 corresponds to L c = Lf/10.) For the 
solutions labeled A, B, and C, Eq is chosen so that U t /UB 
equals 2.4 x 10” 1 , 1.4 x 10 -2 , and 2.2 x 10~ 3 , respectively, 
where 

v 2 

U B = y (10) 

is the energy per unit mass of the background magnetic field. For 
solution A, F) possesses a power-law scaling at ko k <-£ 10 3 ko 

in which k : F k oc k :, ' /2 , as in weak fast- wave turbulence in 
homogeneous plasmas (Cho & Lazarian 2002; Chandran 2005). 
In solution B, k 2 F k is slightly steeper than k~ 2 ' 2 at k ~ 5ko 
but flattens to become oc k -3 / 2 at larger k. On the other hand, 
the fast-wave spectrum is much steeper in solution C than in 
solutions A and B. 


4 We note that if the fast-wave cascade were purely local in k, then Fk would 
decrease exponentially with increasing k in the escape-dominated case. 
However, the full weak-turbulence “collision integral" on the right-hand side 
of Equation (9) includes nonlocal interactions. Thus, the cascade time at large 
k in the escape-dominated limit becomes independent of k and has a value that 
is controlled by the amplitude of the power spectrum near kg. When both the 
cascade time and the escape time are independent of k, one recovers a power 
law whose slope depends on the ratio of these two timescales. 


The differing role of wave escape in these three numerical 
solutions can be understood by comparing the energy cascade 
time r c at k = ko with the escape time r esc . The cascade 
timescale r c (k) can be thought of as the time required for fast- 
wave energy to cascade from wavenumber k to wavenumber 
~3 k. To calculate r c (ko), we carry out a series of numerical 
integrations of Equation (5) with S k — »■ 0 and r esc — » 0. At 
t = 0, we set 


F k 


4U tQ /_*_ 
37 r 3 / 2 k 2 V^o 



( 11 ) 


where U t o is a constant that we vary from one numerical 
calculation to the next. It follows from Equations (4) and (11) 
that U t o is the fast-wave energy per unit mass at t = 0. After 
integrating Equation (5) forward in time, we record the time T c 
at which half of the initial wave energy has been dissipated. The 
results are shown in Figure 2. The solid line in Figure 2 is a fit 
to our numerical results, given by 


T c = 


L4Ma {t£) 


( 12 ) 


This expression can be compared to the energy cascade 
timescale in the inertial range of homogeneous, weak, fast- wave 
turbulence (Chandran 2005), 


T cW 


sin 2 (0)ki>A 




(13) 


where 

Sv k = \J \nk 2 F k (14) 

is approximately the rms amplitude of fast-wave velocity fluctu- 
ations at scale . The inertial range is the range of wavenum- 

bers ko k <5C ko, and kj is the wavenumber at which 
dissipation (in our case hyperviscosity) becomes important. Be- 
cause most of the wave energy is concentrated near k = ko, 

U t ~ (<5u, 0 ) 2 . (15) 

Thus, the value of T c in Equation (12) is comparable to the value 
of r c(k) in Equation (13) when k is set equal to ko, except that the 
sin 2 0 term drops out because 7/ is an effective decay timescale 
integrating over all values of 6. 

In weak fast-wave turbulence in the absence of wave escape, 
k 2 F k oc k~ 3 / 2 , Sv k oc k -1//4 , and r c (k) oc k -1 / 2 (Cho & Lazarian 
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Figure 3. Each symbol corresponds to a numerical solution of Equation (5). The 
dashed lines correspond to the following values of T: 2.16 x 10~ 2 , 4.26 x 10~ 2 , 
1.12 x 10 _1 , and 1.66 x 10 -1 (from lower left to upper right). The resulting 
values of Fdiss/Fo are given in the legend. See the text for details. 


Figure 4. Each symbol corresponds to a steady-state numerical solution of 
Equation (5) with one of the three values of Lf/L c shown in the legend. See the 
text for details. 


2002; Chandran 2005). Because T c (k) oc A -1 / 2 , much of the time 
T c required for energy to cascade from ko to Ad is spent cascading 
from k 0 to ~3&o. Thus, the decay timescale T c is comparable to 
the energy cascade timescale at k = ko', i.e., 

Tc(ko) ~ T c . (16) 


On physical grounds, the energy cascade time at k = ko depends 
primarily on the values of ko, £/ t , and i; A (which determine the 
strength of nonlinear wave-wave interactions at A = ko) and not 
on whether the turbulence is decaying or forced. We thus take 
Equation (16) to apply to the forced-turbulence calculations that 
we present below, but we replace U t o in Equation (12) with U t , 
the instantaneous value of the fast-wave energy per unit mass. 
Equations (12) and (16) then allow us to write 



Tc(ko) 


(17) 


where 


LfUt 
L c u B ■ 


(18) 


Here, we have re-expressed ko in terms of the correlation length 
of the fast-wave fluctuations, which we take to be (see the 
Appendix) 



(19) 


We note that the quantities on the right-hand side of 
Equation (18) are macroscopic parameters that are potentially 
measurable from observations or MHD simulations of flares. 

Returning to the numerical results in Figure 1, when 
7'esc/t'c(ko) is small (as it is for solution C), waves at ko escape be- 
fore their energy cascades to larger k. This causes A) to steepen 
and reduces the fraction of the fast-wave energy that is dissi- 
pated at large k (i.e., transferred to electrons via TTD in a flare). 
On the other hand, when the ratio in Equation (17) is >1 (as it 
is for Solution A), fast-wave energy cascades from ko to larger 
wavenumbers before the waves escape, and /■)■ approaches the 
scaling found in homogeneous turbulence simulations without 
wave escape. 

In Figure 3, we summarize results from 1 8 different numerical 
solutions to Equation (5). To obtain these solutions, we start by 
fixing the value of Ef/E c at 5, 10, 20, 30, or 40. We then vary 
Eq until T has one of four values: 2.16 x 10 -2 , 4.26 x 10~ 2 , 
1.12 x 10~', or 1.66 x 10 _1 . In the figure, we plot one symbol 
for each of the solutions, with a different symbol for each of the 


four values of T. After each solution reaches steady state, we 
calculate the dissipation power 

Ediss = J sin 2 (d)k s vFi < d^k (20) 

in that solution. Since the dissipation in our simulation is a 
proxy for the transfer of large- A' wave energy to electrons via 
TTD, the ratio Ediss/ Eq is a measure of the efficiency of electron 
acceleration. As shown in the legend, the value of this ratio is 
roughly the same for all simulations with the same value of 
T , i.e., along each of the four curves. This demonstrates that 
Ediss/Eo depends almost solely on T, i.e., on r esc /T c (Ao). 

In Figure 4, we plot Edi ss /£o for values of T ranging from 
8.6 x 10” 3 to 5.0. For these calculations, we employed three 
different values for the quantity Ef/L c : 5, 10, and 20. As the 
figure shows, all three sets of calculations produce essentially 
the same plot, again supporting our assertion that Ediss/Eo 
in our model is determined almost exclusively by the single 
quantity T. The solid line in this plot corresponds to the function 


Ediss 



0 . 012 ^ 


-20+0.02T- 1 


( 21 ) 


and provides a reasonable fit to our numerical results. 

For fast waves propagating in low-/! plasmas, half of the wave 
energy is in the velocity fluctuations and half is in the magnetic- 
field fluctuations. The total fast-wave energy is thus twice the 
energy in the fluctuating magnetic field. If SB is defined to 
be the rms amplitude of the magnetic-field fluctuations, then 
Ui/Ub — 2(S B / Bo) 2 . We can thus rewrite Equation (18) as 


T = 




( 22 ) 


In Figure 5, we use Equations (21) and (22) to plot £diss/£o 
as a function of SB /Bo for three choices of the ratio Ef/L c . 
The results show that the rate of energy dissipation (i.e., of 
energy transfer to electrons via TTD) falls off rapidly below 
rms fast-wave amplitudes SB /Bo of about 5%-10%, depending 
on the size of the flare region relative to the correlation length 
(outer scale) of the turbulence. In the strong-turbulence limit 
of larger fast-wave amplitudes, on the other hand, the energy 
dissipated at shorter wavelengths (transferred to electrons via 
TTD) approaches 100% of the input wave energy. 
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Figure 5. Each of the three curves is a plot of Equation (21) for a different 
choice of Lf/L c . See the text for details. 


4. CONCLUSION 

We have developed a leaky-box model for fast-wave turbu- 
lence in solar flares that accounts for wave escape from the solar- 
flare acceleration region in an approximate way. In this model, 
we determine the fast-wave power spectrum Fk by solving the 
wave kinetic equation for fast waves from weak turbulence the- 
ory (Chandran 2005, 2008) after modifying this equation in the 
following ways. First, we set the amplitudes of the other wave 
types (the slow magnetosonic wave and Alfven wave) to zero, 
thereby neglecting interactions between fast waves and other 
wave types. Second, we add a homogeneous loss term to account 
for wave escape. Third, we add a source term S(k ) to model the 
generation of fast waves by the disordered flows that are pre- 
sumed to arise when reconnection outflows encounter closed 
magnetic loops lower in the solar atmosphere. This source term 
is isotropic in wavenumber space, peaks at a wavenumber ko, 
and injects fast-wave energy per unit mass at the rate Eq, which 
is an adjustable parameter. Fourth, we add a hyperviscous dis- 
sipation term, so that virtually all the energy that cascades to 
large k is dissipated. The energy dissipation rate per unit mass 
£diss is a proxy for the power (per unit mass) that is available to 
accelerate particles through wave-particle interactions at large 
wavenumbers. 

We numerically integrate this modified wave kinetic equation 
forward in time until a steady state is reached, varying Eq and 
the quantity koLf, where Lf is the size of the flare acceleration 
region. We find that the acceleration efficiency, E&ss/Eq, de- 
pends on a single quantity, T (defined in Equation (18)), which 
is approximately the ratio of the escape time r esc to the energy 
cascade time at ko, which is denoted r c (ko). When r esc <SC T c (ko), 
almost all of the fast-wave power that is injected at wavenum- 
ber ko escapes the flare acceleration region before cascading to 
small scales. On the other hand, if T esc > t c (ko), then a sizable 
fraction of the fast-wave power that is injected at ko cascades 
to large k and dissipates before leaving the flare acceleration 
region. 

To gain further insight into the effects of wave escape on 
electron acceleration by fast waves in solar flares, it will be 
important to obtain stronger constraints on if/ Uq and Lf/ L c 
in flares. One way that such constraints could be obtained 
is through numerical simulations of magnetic reconnection, 
outflows, and large-scale magnetic structure in solar flares. 
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APPENDIX 

RELATION BETWEEN WAVENUMBER AND 
CORRELATION LENGTH 

We consider a superposition of waves, all with the same 
wavenumber k, but with a distribution of wavevector directions 
that becomes isotropic in the limit that the number of waves is 
increased toward infinity. We take the phase of each wave to have 
a random additive phase constant that is uniformly distributed 
between 0 and 2jt. If fix) is some quantity that is modulated by 
waves — e.g., a component of the magnetic field — then we can 
write 

N 2 N 

f(x) = sin(*, 7 • x + i/fij), (Al) 

i=o j=i 

where k t j — £[(1 — /xj) l/2 coscpjx + (1 — /r 2 ) 1 / 2 sin cpjy + fijz]. 
Hi = 2 i/N — 1, <pj = 2 itj/N, and i Jr^ is the random, additive 
phase constant. We set 


A (k ,j ) — A ; v , 


(A2) 


where is independent of kjj but decreases as N increases 
so that (|/(x)| 2 ) remains finite as N — ► oo. Because A ^ is 
independent of k,, , the power spectrum of / becomes isotropic 
as N -» oo. 

The autocorrelation function of f(x) is 


( f(x 0 )f(XQ + X )) 
(/(*o)/(*o)> 


(A3) 


where (. . .) denotes an ensemble average — i.e., an average over 
the random phases of the waves. Making use of the relations 


(sin(i/ ly ) sin(i/f m „)) = (cos(i/r, 7 ) cos(i/f m „)) = 




jn 


(A4) 


and 

(sin((/r, 7 ) cos(i/f m „)) = 0, (A5) 


we find that 


a 2 n IN 

a n \ " \ " 


o 

K 

o 

+ *)} = X! cos(k ij ■ 

<=0 7=1 

(A6) 

and 

</(*o)/(*o)) = X, 

(A7) 

where 

X = lim N(N+1)A 2 n . 

N^-oo 

(A8) 

In the limit N - 

-> oo, Cix) depends only on \x 

. We thus 

define 

C(x ) = lim C(x). 

N—>oo 

(A9) 


Without loss of generality, we set x = xz to obtain 


This work benefited from valuable discussions with our 
colleagues in a NASA Living-With-a-Star Focused-Science- 
Topic team working on “Flare Particle Acceleration Near 


/ f 1 y sin(Ct) 

hm {f{x 0 )f(x 0 +x)) = - / dixcos(Hkx) = 

jv-»oo 2 J_i kx 

(A10) 
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Equations (A7) and (A10) imply that 


PONGKITIWANICHAKUL ET AL. 


sin(Cx) 

kx 


(All) 


We define the correlation length of f(x), L c , through the relation 


C(L c ) = e~ 1 . 


This gives 


2.2 



(A 12) 
(A 13) 
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